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We study the magnetic excitation spectra in the paramagnetic state of BaFe2As2 from the ab 
initio perspective. The one-particle excitation spectrum is determined within the combination of the 
density functional theory and the dynamical mean- field theory method. The two-particle response 
function is extracted from the local two-particle vertex function, also computed by the dynamical 
mean field theory, and the polarization function. This method reproduces all the experimentally 
observed features in inelastic neutron scattering (INS), and relates them to both the one particle 
excitations and the collective modes. The magnetic excitation dispersion is well accounted for by 
our theoretical calculation in the paramagnetic state without any broken symmetry, hence nematic 
order is not needed to explain the INS experimental data. 



Neutron scattering experiments provide strong con- 
straints on the theory of iron pnictides. Both the lo- 
calized picture and the itinerant picture of the mag- 
netic response have had some successes in accounting or 
even predicting aspects of the experiments. Calculations 
based on a spin model with frustrated exchange con- 
stants [HQ or with biquadratic interactions H described 
well the neutron scattering experiments 0, HI] . The itin- 
erant magnetic model, based on an random phase ap- 
proximation (RPA) form of the magnetic response, uses 
polarization functions extracted from density functional 
theory (DFT) [6] or tight binding fits 0^ and produces 
equally good descriptions of the experimental data. 

Furthermore, DFT calculations predicted the stripe 
nature of the ordering pattern [To| and the anisotropic 
values of the exchange constants which fit well the spin 



wave dispersion in the magnetic phase The tight 

binding calculations based on DFT bands also predicted 
the existence of a resonance mode in the superconducting 
state 



121. 



In spite of these successes, both itinerant and localized 
models require significant extensions to fully describe the 
experimental results. DFT fails to predict the observed 
ordered moment [ll|. Furthermore, adjusting parame- 
ters such as the arsenic height to reproduce the ordered 
moment, leads to a peak in the density of states at the 
Fermi level instead of the pseudogap, which is ob- 
served experimentally. The localized picture cannot de- 
scribe the magnetic order in the FeTe material without 
introducing additional longer range exchange constants. 
Given that this material is more localized than the 122, 
the exchange constants would be expected to be shorter 
range. Furthermore, fits of the INS data require the use 
of anisotropic exchange constants well above the mag- 
netic ordering temperature [13| . However no clear phase 
transition to a nematic phase in this range has been de- 
tected. 

In this Letter, we argue that the combination of den- 
sity functional theory and dynamical mean field theory 
(DFT+DMFT) provides a natural way to improve both 
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Figure 1: (Color online) The Feynman diagrams for the 
Bethe-Salpeter equation. It relates the two-particle Green's 
function (%) with the polarization (x°) and the local irre- 
ducible vertex function (T\ r r c ). The non-local two-particle 
Green's function is obtained by replacing the local propagator 
by the non-local propagator. 



the localized and the itinerant picture, and connects the 
neutron response to structural material specific informa- 
tion and to the results of other spectroscopies. 

We compute the one-particle Green's function using 
the charge self-consistent full potential DFT+DMFT 
method, as implemented in Ref. Q, based on Wien2k 
code [15]. We used the continuous-time quantum Monte 
Carlo (CTQMC) [l6|,[l7| as the quantum impurity solver, 
and the Coulomb interaction matrix as determined in 
Ref. [H. The dynamical magnetic susceptibility x(q,a;) 
is computed from the ab initio perspective by extract- 
ing the two-particle vertex functions of DFT+DMFT 
solution T\ r r c [19]. The polarization bubble \° is com- 
puted from the fully interacting one particle Greens func- 
tion. The full susceptibility is computed from x° an d the 



two-particle irreducible vertex function T\ r J c , which is as- 
sumed to be local in the same basis in which the DMFT 
self-energy is local, implemented here by the projector to 
the muffin-thin sphere 



14j . In order to extract T^, we 



employ the Bethe-Salpeter equation (see Fig.[T]) which re- 
lates the local two-particle Green's function (xioc), sam- 
pled by CTQMC, with both the local polarization func- 

tion (X? OC ) and F 
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depends on three Matsubara frequencies (iv, iv '; iw), 
and both the spin (01-4) and the orbital (oti-4) indices, 
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Figure 2: (Color online) (a) The constant energy plot 
of the theoretical dynamical structure factor, S(q, ou) 

(=- x _i2 , / J k B r ) a ^ different energies (50meV, 75meV, 
125meV, and 150meV) in the paramagnetic state (T=386 K) 
of BaFe2As2 as a function of momentum q =(H,K,L). L is 
here fixed at 1. (b) The corresponding inelastic neutron scat- 
tering data from Ref. [l3|. 



which run over 3d states on the iron atom. T is the 
temperature. 

Once the irreducible vertex V\ r r c is obtained, the mo- 
mentum dependent two-particle Green's function is con- 
structed again using the Bethe-Salpeter equation (Fig.[TJ 
by replacing the local polarization function x? oc by the 
non-local one X q ,iu; : 
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Finally, the dynamic magnetic susceptibility x(q 7 ^j) is 
obtained by closing the two particle green's function with 
the magnetic moment /i = /is(L + 2S) vertex, and sum- 
ming over frequencies (iv,iv'), orbitals (ai-4), and spins 
(cri-4) on the four external legs 
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The resulting dynamical magnetic susceptibility is 
obtained in Matsubara frequency (ioj) space and it 
needs to be analytically continued to real frequencies 
(x(q, oj)). For the low frequency region, on which we 
concentrate here, the vertex Y\ r r c is analytically con- 
tinued by a quasiparticle-like approximation. We re- 
place the frequency dependent vertex with a constant, 
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Figure 3: (Color online) S(q, oj) along the special path in 
Brillouin zone marked by red arrow in the inset on the right. 
The inset shows the body-centered tetragonal (black line) and 
the unfolded (blue line) Brillouin zone. Black dots with error 
bars correspond to INS data from Ref. 13. The white dashed 
line shows the isotropic Heisenberg spin wave dispersion. 
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x(<7, ioj = 0) = xfe^ 7 = 0). This vertex U however 
retains important spin and orbital dependence. 

Fig. [2ja) shows the calculated constant energy plot of 



the dynamical structure factor, S(q, oj) in the paramag- 
netic state of BaFe2As2- Our theoretical results are cal- 
culated in the unfolded Brillouin zone of one Fe atom per 
unit cell, because magnetic excitations are concentrated 
primarily on Fe atoms, therefore folding, which occurs 
due to the two inequivalent arsenic atoms in the unit cell, 
is not noticeable in magnetic response [7]. For compari- 
son we also reproduce in Fig. EJb) the INS experimental 
data from Ref. [l3|. At low energy (around o;=50meV), 
the theoretical S(q, oj) is strongly peaked at the order- 
ing wave vector (H,K,L) = (1, 0, 1) and it forms a clear 
elliptical shapes elongated in K direction. The elonga- 
tion of the ellipse increases with energy (cj=75meV) and 
around oj =125meV the ellipse splits into two peaks, one 
peak centered at (1,0.4, 1) and the other at (1, —0.4, 1). 
At even higher energy (oj ~150meV) the magnetic spec- 
tra broadens and peaks from four equivalent wave vec- 
tors merge into a circular shape centered at wave vector 
(1,1,1). At even higher energy (230meV, not shown in 
>^the figure) the spectra broadens further, and the peak 
becomes centered at the point (1,1,1). These trends are 
all in good quantitative agreement with INS data from 

Fig. EDO . 

In Fig. [3j we display a contour plot of the theoret- 
ical S(q, oj) as a function of frequency uo and momen- 
tum q along the special path in the unfolded Brillouin 
zone, sketched by a red line in the right figure. At low 
energies (uo <80meV), S(q, oj) is mostly concentrated in 
the region near the ordering vector (1,0,1). Consistent 
with the elongation of the ellipse along the K direction in 
Fig. [2j the low energy (oj <80meV) bright spot in Fig. [3] 
is extended further towards (1,1,1) direction but quite 
abruptly decreases in the (0,0,1) direction. The mag- 
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netic spectra in the two directions (1,0, 1) —> (0,0,1) and 
(1,0,1) — » (1,1,1) are clearly different even at higher 
energy uj > lOOmeV. The peak position is moving to 
higher energy along both paths, but it fades away very 
quickly along the first path, such that the signal prac- 
tically disappears at (0.5,0,1). Along the second path 
(1,0,1) — »• (1,1,1), there remains a well defined ex- 
citation peak for which the energy is increasing, and 
at (1,1,1) reaches the maximum value of ~ 230meV. 
Continuing the path from (1,1,1) towards (0,0,1) the 
peak energy decreases again and it fades away around 
(0.5, 0.5, 1). The black dots display INS data with errors 
bars from Ref. 13. Notice a very good agreement between 
theory and experiment. 

The white dashed line in Fig.[3]represents the spin wave 
dispersion obtained for the isotropic Heisenberg model 
using nearest neighbor J\ and next nearest neighbor J2 
exchange constants and performing the best fit to INS 
data. This fit was performed in Ref. [l3|. The magnetic 
excitation spectra of an isotropic Heisenberg model show 
a local minimum at the wave vector q = (1,1,1), which 
is inconsistent with our theory and with the experiment. 
To better fit the experimental data with a Heisenberg- 
like model, very anisotropic exchange constants need to 
be assumed [13|, which raised speculations about possi- 
ble existence of nematic phase well above the structural 
transition of BaFe 2 As 2 . Since the DFT+DMFT results 
can account for all the features of the measured mag- 
netic spectra without invoking any rotationally symme- 
try breaking the presence of nematicity in the param- 
agnetic tetragonal state at high temperature is unlikely. 

In Fig. Ufa) we show constant frequency cuts in the 
K direction (from (1,-1,1) through (1,0,0) to (1,1,1)) 
of S(q, uj) displayed in Fig. [3l For comparison we also 
show the corresponding INS measurements from Ref. [HI 
as red circles in Fig. |4jb) and (c). At cj=20meV, the 
spectrum has a sharp peak centered at the ordering vec- 
tor (1,0, 1). At cj=50meV, the spectrum still displays a 
peak at (1, 0, 1) but the intensity is significantly reduced. 
With increasing frequency gj, the peak position in S(q, uj) 
moves in the direction of (1,1,1), and at 128meV peaks 
around (1,0.4, 1). The shift of the peak is accompanied 
with substantial reduction of intensity at ordering wave 
vector (1, 0, 1). At even higher energy of 250meV only a 
very weak peak remains, and it is centered at the wave 
vector (1,1,1). The position of peaks as well as their fre- 
quency dependence is in very good agreement with INS 
experiments of Ref. LL3| displayed in Fig. HJ^b) and (c). 

Fig. [5fa) resolves the dynamical magnetic suscep- 
tibility x °f Eq. [3] in the orbital space Xa — 

fi z aia3 fi z a2(J4 X- 1,^2(^,^0 such 
0-30-4 CK0-3, po-4 

that x — ^2 a Xa- At the magnetic ordering vector 
q=(l,0, 1), Xa increases sharply with frequency near 
uj = for all orbitals and is strongly suppressed above 
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Figure 4: (Color online) (a) The wave vector K dependence 
(H=1,L=1) of S(q, uj) at several frequencies, (b) The corre- 
sponding INS data at cj=19meV and (c) 128meV reproduced 
from Ref. [l3|. The red circles correspond to the paramagnetic 
state at T=150K and the blue diamonds to the magnetic state 
at T=7K. 



lOOmeV reaching the maximum around 20meV. At this 
wave vector, the dominant contributions at low energy 
come from the d xy and the d yz orbitals. The magnetic 
susceptibility at q=(0, 1,1) in Fig. [5fa) shows the same 
trend as orbit ally resolved spectra at q=(l,0, 1), except 
that d xz and d yz switch their roles due to the C4 sym- 
metry of the Fe square lattice. 

These dominant orbital contributions to x are a l so 
reasonably captured in the polarization bubble x° ( n °t 
shown here), hence these excitations could be understood 
in terms of the Fermi surface nesting. The orbital re- 
solved Fermi surface is displayed in Fig. Efb) at both the 
T-plane and the Z-plane. Most of the weight in x° comes 
from the diagonal terms, i.e., Xa hence the Fermi sur- 
faces with the same color in Fig. [5fb) but separated by 
the wave vector (1, 0, 1) give dominant contribution. The 
intra-orbital d yz low energy spectra comes mostly from 
the transitions between the green parts of the hole pocket 
at T and the green parts of the electron pocket at A, 
marked with green squares (□) in Fig. [51(b). Since the 
electron pocket at A is elongated in H direction, the nest- 
ing condition occurs mostly in the perpendicular K di- 
rection, hence the elliptical excitations at low energy in 
Fig. [2] are elongated in K but not in H direction. The 
intra-orbital d xy transitions are pronounced between the 
electron pocket at M' and the hole pocket at R, as well 
as between the electron pocket at A! and the hole pocket 
at X (marked with red Q)- This large spin response at 
(1, 0, 1) gives rise to the low energy peak in Fig. [3] 

We note that the particle-hole response, encoded in 
polarization bubble x°? is especially large when nesting 



(a) 



3 



A (H,K,D = (1,0,1) 


A (H,K,L)=(0,1,1) 






(H,K,L)=(0,0,1) 


— z 2 

— x 2 -y 2 

— xz 

— yz 

— xy 


(H,K,D = (1,1,D 









100 200 300 400 100 200 300 400 500 




10 1 



H(r.l.u.) H(r.l.u.) 




Figure 5: (Color online) (a) The Fe d orbital resolved dy- 
namical magnetic susceptibility at T=386K for distinct wave 
vectors q=(l,0,l), (0,1,1), (0,0,1), and (1,1,1). Different 
colors correspond to different orbital contributions, (b) The 
Fermi surface in the Y and Z plane at T=73K colored by 
the orbital characters: cfcz(blue), d yz (green), and d xy (ied). 
The small symbols mark the regions in the Fermi surface, 
where nesting for the wave vector q=(l, 0, 1) is good, (c) The 
zoom-in of A(k,oj) along the path marked by black dashed 
line in Fig. [5jb). The green open circles indicate the two rel- 
evant bands of d xy character which give rise to the peak in 
magnetic excitation spectra near 230meV at q=(l, 1,1). (d) 
A(k, uo) computed in one Fe atom per unit cell at T=73K. The 
green arrows mark the same bands which give rise to 230meV 
peak. The DFT bands are overlayed by white dashed lines. 
The blue arrows mark corresponding DFT bands of d xy char- 
acter. 



occurs between an electron pockets and a hole pocket, 
because the nesting condition extends to the finite fre- 
quency, and is not cut-off by the Fermi functions. 

The low energy magnetic excitations at wave vectors 
q = (0,0,1) and q = (1,1,1) can come only from 
electron-electron or hole-hole transitions, hence both re- 
sponses are quite small, as seen in Fig. [5fa). While the 
magnetic response at q = (0, 0, 1) is small but finite, the 
spin response at q = (1,1,1) is almost gapped. This is 
because the hole-hole transitions from V to R or electron- 
electron transitions from M to A' do not involve any 
intra-orbital transitions, and hence are even smaller than 
transitions at the wave vector (0, 0, 1). 

At finite energy transfer, the spin excitations come 
from electronic states away from the Fermi energy, and 
can not be easily identified in the Fermi surface plot. 
Hence it is more intriguing to find the dominant contri- 
bution to the peak at uo « 230meV and q = (1,1,1). This 
peak gives rise to the 230meV excitations at (1,1,1) in 
Fig. 03 A large contribution to this finite frequency exci- 
tation comes from a region near the two electron pockets 
at M and A' marked with black dashed line in Fig. 03(b). 
We display in Fig. [5](c) the one electron spectral function 
across these dashed lines in the Brillouin zone to show 
an important particle hole transition from the electrons 
above Fermi level at the M point and the flat band at 
-200meV around the A' point, both of d xy character. We 
note that due to large off diagonal terms in the two parti- 
cle vertex T, all orbital contributions to x develop a peak 
at the same energy, although only d xy orbital displays a 
pronounced peak in x°- 

Fig. 03(d) displays the one electron spectral function 
in a path through the Brillouin zone, corresponding to 
one Fe atom per unit cell. Within DFT+DMFT the 
quasi-particle bands are renormalized by a factor of 2- 
3 compared to the corresponding DFT bands (white 
dashed lines). The green arrow marks the d xy band 
which contributes to the peak in 5(q, uo) near 230meV 
and q = (1,1,1). In DFT calculation, this d xy intra- 
orbital transition is also present, but occurs at much 
higher energy of the order of 400-600meV, marked by 
blue arrows. The over-estimation of the peak energy at 
q=(l, 1,1) was reported in LSD A calculation of Ref. 0. 

In this Letter, we have extended the DFT+DMFT 
methodology to compute the two particle responses in 
a realistic multi-orbital DFT+DMFT setting. With the 
same parameters which were used to successfully describe 
the optical spectra and the magnetic moments of this 
material [20], we obtained a coherent description of the 
experimental neutron scattering results. Our theory ties 
the magnetic response to the fermioloy of the model, and 
quantifies the departure from both purely itinerant and 
localized pictures. 
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